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We consider a simple configuration for realizing one-dimensional slow-light metamaterials with 
large bandwidth-delay products using stub-shaped Fabry-Perot resonators as building blocks. Each 
metaatom gives rise to large group indices due to a classical analog of the dressed-state pic¬ 
ture of electromagnetically-induced transparency. By connecting up to eight metaatoms, we hnd 
bandwidth-delay products over unity and group indices approaching 100. Our approach is quite 
general and can be applied to any type of Fabry-Perot resonators and tuned to different operating 
wavelengths. 


Electromagnetically-induced transparency (EIT) oc¬ 
curs iu atomic systems due to the destructive quau- 
tum iuterfereuce of coupled atomic trausitiousii^. The 
EIT pheuomeuou readers au otherwise optically opaque 
medium trauspareut iu a uarrow spectral regiou, which, 
due to the Kramers-Krouig relatious, is accompauied by 
stroug material dispersiou, resultiug iu a siguificaut re- 
ductiou of the speed of light^ or eveu stoppiug light^i^. 
For all-optical processiug systems, coutrolliug the speed 
of light is au esseutial ability, especially for desiguiug op¬ 
tical storage systems, such as optical buffers^. Uufortu- 
uately, atomic EIT ofteu requires extreme euviroumeu- 
tal couditious, such as cryogeuic temperatures, which are 
uot suitable for chip-based systems. Iu coutrast, classical 
aualogs of EIT iu metal-based plasmouic structures^-— 
pave the way for exploitiug the reductiou of speed of light 
ou the uauoscale. Such extreme coutrol of light is possi¬ 
ble due to the coupliug of light to surface plasmous prop- 
agatiug at metal-dielectric iuterfaces, thereby squeeziug 
the light beyoud the diffractiou limifii^. 

Iu the quautum mechauical descriptiou of atomic EIT, 
two theoretical iuterpretatious based ou either the bare 
states or the dressed states are available. While both 
are equivaleut iu the quautum descriptiou, their classi¬ 
cal aualog differ siguificautly. Iu particular, the classi¬ 
cal EIT aualog of the bare states iuterpretatiou requires 
the use of radiaut aud subradiaut plasmouic modesi"— 
while the dressed states picture is based ou detuued radi¬ 
aut resouator s^^d^d^ . The dressed-states aualog has re- 
ceutly received substautial atteutiou, especially by utiliz- 
iug detuued resouators coupled to metal-iusulator-metal 
(MIM) waveguidesi^ii^-— . The gap surface plasmou 
mode supported by the MIM waveguide offers a good 
compromise betweeu mode coufiuemeut aud propagatiou 
leugth, makiug it au ideal caudidate for slow-light de¬ 
vices. 

Oue of the parameters for characteriziug a slow-light 
device is the miuimal group velocity or, equivaleutly, 
the largest group iudex Ug = c/v^, which cau be attaiued. 
However, a large group iudex is ouly practical if avail¬ 
able iu a wide frequeucy baudwidth A/. Heuce, au im- 
portaut figure-of-merit for slow-light applicatious is the 
baudwidth-delay product r A/, where r = v^l is the delay 


experieuced by au optical pulse over a propagatiou leugth 
To eusure the delay of a pulse of duratiou ^ A/“^, 
we must have rAf > I. Iu this work, we preseut a sim¬ 
ple coufiguratiou for geueratiug over-uuity baudwidth- 
delay products by creatiug oue-dimeusioual metamate¬ 
rials usiug two detuued Fabry-Perot resouators as our 
metaatom. We study detuued stub resouators aperture- 
coupled to a MIM waveguide to realize the dressed-state 
aualog of atomic EIT with accompauyiug large group iu¬ 
dex. By joiuiug up to eight metaatoms, we show the ex¬ 
pected couvergeuce of both the real aud imagiuary parts 
of the metamaterial effective refractive iudex. We fiud a 
group iudex approachiug 100 aud baudwidth-delay prod¬ 
uct of 1.2 with a reasouable trausmittauce of 10%. 

We cousider a silver-glass-silver waveguide, where the 
refractive iudex of glass is 1.45, while the dielectric fuuc- 
tiou for silver is iuterpolated from tabulated values2^. 
The width of the waveguide is fixed to re = 100 um 
throughout this study. Usiug a port excitatiou scheme 
iu a commercial fiuite-elemeut frequeucy-domaiu solver 
(COMSOL Multiphysics 5.0), we lauuch the gap sur¬ 
face plasmou mode of the MIM waveguide aud detect 
the trausmittauce at au output port positioued after the 
stub resouators. Before discussiug the system cousistiug 
of two detuued stub resouators, we first optimize the case 
of a siugle stub resouator of leugth Li, which is coupled 
through au aperture of width d aud gap g to the MIM 
waveguide, as illustrated iu Fig. [D^a). Iu Fig. [TJb) we 
vary the leugth of the stub for fixed aperture dimeusious 
of ^ = d = 50 um. For Li < 500 um, we uotice a siugle 
dip iu the trausmittauce spectra due to the excitatiou of 
au auti-symmetric Fabry-Perot mode supported by the 
stub, which exhibits a staudiug wave patteru with two 
autiuodes iu the maguetic field (uot depicted). The traus¬ 
mittauce dip is due to the destructive iuterfereuce of the 
bus MIM wave aud the reflected wave from the stub res¬ 
ouator. For Li> 500 um au additioual trausmittauce dip 
at shorter waveleugth shows up siuce the stub resouator 
uow also supports a symmetric Fabry-Perot mode with 
three autiuodes. For iucreasiug stub leugth, both Fabry- 
Perot modes redshift aud broadeu. By varyiug the stub 
leugth, we cau tuue the resouauce of a giveu Fabry-Perot 
mode to a desired waveleugth. Iu the subsequeut aualy- 
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FIG. 1. (a) Schematic illustration of silver-glass-silver waveg¬ 
uide aperture-coupled (with width d and gap g) to a single 
stub-shaped Fabry-Perot resonator of length Li. Transmit¬ 
tance through waveguide (color plot) as a function of wave¬ 
length and (b) stub length Li, (c) aperture width d, and (d) 
aperture gap g. The following parameters are used, unless 
they are varied in the plot: ic = 100 nm, d = ^ = 50 nm, and 
Li = 500 nm. 


sis, we choose Li = 500 nm since both modes are then 
present in the considered wavelength range. In Fig. [T]^c) 
and (d) we study the effect of the aperture width d and 
gap g, respectively. As an important feature, we see that 
decreasing the aperture width provides a narrower trans¬ 
mittance dip since the reflection from the stub surface 
next to the aperture increases. Put in other words, the fi¬ 
nesse of the stub-shaped Fabry-Perot resonator increases 
with decreasing aperture width. As shown in earlier stud¬ 
ied, a spectrally narrow transmittance dip is beneficial 
since it produces a strong EIT effect with resulting large 
group index. However, we find that for too small val¬ 
ues of d the transmittance minima begin to increase (i.e., 
the transmittance does not reach zero at the resonance 
wavelengths) due to the decrease in coupling between the 
MIM waveguide and the stub resonator. Hence, a trade¬ 
off between the transmittance minimum and the spectral 
width of the transmittance dip (i.e., quality factor) ex¬ 
ists, which can be controlled through the aperture width. 
Finally, as shown in Fig.lTJd) increasing the aperture gap 
(for a fixed aperture width) only produces a redshift of 
the transmittance dip (and no significant change in the 
spectral width) as a consequence of additional phase ac¬ 
cumulation in the aperture, which slightly modifies the 
resonance condition of the aperture-coupled stub system. 
Based on this analysis, we keep the same value for the 


aperture gap (i.e., = 50 nm) but decrease the width 

to d = 30 nm for a narrow transmittance dip which still 
reaches zero at the Fabry-Perot resonance wavelengths. 

Next, we add an additional stub of length L 2 on the 
opposite side of the waveguide, as illustrated in Fig.[2]^a). 
The two apertures, which couple to each of the stub 
resonators, have the same dimensions. We will later 
consider the double-stub system as our metaatom (with 
unit size a) when building our one-dimensional slow-light 
metamaterial. For now, we study a single cell. To gener¬ 
ate the EIT effect based on the dressed state interpreta¬ 
tion we need two detuned resonator s^^b^ ^ which can be 
achieved by simply varying the length of the second stub 
resonator, while keeping the widths of the two stubs iden¬ 
tical (the opposite approach will also give rise to detuned 
resonators). Eigure[2fb) shows transmittance spectra of 
the double-stub system for L 2 values which are slightly 
larger than Li {Li = 500 nm). As the difference be¬ 
tween the lengths of the two stubs L 2 — Li increases, 
the detuning of the two resonators increases, creating an 
EIT window in the spectral regions of the transmittance 
dips. We note that the long-wavelength transmittance 
dip requires a larger detuning to create the EIT window, 
since the dip is spectrally wider. In fact, to have a suf¬ 
ficiently large transmittance maxima in the EIT window 
(^ 80%) we find that L 2 should be at least 550 nm. In 
Fig. HJc-e) we show the magnetic field profile of the EIT 
phenomenon for three different wavelengths for the case 
of L 2 = 550 nm. Eor wavelengths shorter [Eig. [2fc)] and 
longer [Eig. EJe)] than the EIT window, the two stub 
resonators are excited independently, producing in both 
cases a weak transmission through the waveguide. At the 
transmittance maximum in the EIT window [Eig. Efd)], 
both stub resonators are excited simultaneously but with 
opposite phase, thereby cancelling the destructive inter¬ 
ference normally created by the individual stubs. A sim¬ 
ilar picture is seen for the low-wavelength EIT window 
(with three antinodes in the magnetic field profile of each 
stub instead of two), showing that the EIT effect created 
due to the detuned resonators is independent of the stub 
mode. 

In Eig. [3] we focus on the long-wavelength EIT window 
[Eig. [3](a)] and determine the effective refractive index 
rieff [Eig.[3l^b-c)] along with the group index rig [Eig.[3l^d)]. 
The real part of the effective index is calculated from the 
phase difference of the magnetic fields Acj) at positions be¬ 
fore and after the metaatoms as Re(neff) = Acjr/^Nako)^ 
where N is the number of metaatoms and ko = 2'kcI\ 
is the vacuum wave vector— The imaginary part 
of the effective index is determined from the transmit¬ 
tance T as Im(neff) = — ln(T)/(2A^a/co), assuming an 
exponential decay of the power of the fields. Einally, the 
group index is given as rig = Re(neff) — the 

spectral regions of normal dispersion <0. To 

avoid unphysical large variations in the group index, we 
perform a linear fit of the real part of the effective in¬ 
dex in the regions of normal dispersion [shown as dashed 
lines in Eig. [3][b)]. Eigure [3] shows that a single double- 






























3 




FIG. 2. (a) Schematic illustration of silver-glass-silver waveg¬ 
uide aperture-coupled (with width d and gap g) to two stub¬ 
shaped Fabry-Perot resonators of lengths Li and L 2 . The 
metaatom cell size a is also shown, (b) Transmittance through 
the waveguide (color plot) as a function of wavelength and sec¬ 
ond stub length L 2 . The hxed parameters are: re = 100 nm, 
d = 30 nm, ^ = 50 nm, and Li = 500 nm. (c-e) Out-of- 
plane magnetic held component at the resonance wavelength 
of the short stub, the wavelength of maximum transmittance 
in the FIT-like window, and the resonance wavelength of the 
long stub, respectively. The length of the second stub is 
L 2 = 550 nm. 


stub metaatom produces a group index larger than 60 
in a bandwidth of approximately 50 nm [determined as 
the full-width at half-maximum (FWHM) of the trans¬ 
mittance peak], with a transmittance maximum close to 
80%. 

The metamaterial description is only applicable for 
metaatoms signihcantly smaller than the operating wave¬ 
length, i.e., a ^ A, setting an upper limit to our cell 
size. Simultaneously, too small cell sizes give rise to un¬ 
desired coupling between the metaatoms. Despite these 
two limits the cell size is still quite tunable, and we 
find that a cell size of a = 200 nm provides a good 
trade-off. Figure EJb-c) shows that the effective index 
quickly converges with increasing number of cells, sup¬ 
porting the metamaterial interpretation. Additionally, 
the group index [Fig. [3][d)] also increases with the num¬ 
ber of cells, reaching close to a two-order decrease in 
the speed of light with only 8 cells. However, due to 
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FIG. 3. (a) Transmittance T through one-dimensional slow- 
light metamaterial consisting of one to eight double-stub 
metaatoms, (b-c) Real and imaginary parts of the effective 
index neff of the metamaterial, respectively. In the spectral 
regions of normal dispersion, a linear ht of the real part of 
the effective index is performed (dashed lines), (d) Group 
index rig determined from the linear hts of the real part of 
the effective index. The metaatom consisting of two detuned 
aperture-coupled stubs is sketched in Fig. [2l[a) with the fol¬ 
lowing parameters: w = 100 nm, d = 30 nm, ^ = 50 nm, 
Li — 500 nm, L 2 = 550 nm, and a = 200 nm. 


the increase in propagation length I = Na the FWHM 
and maximum transmittance decrease. As discussed ear¬ 
lier, an important figure-of-merit for any slow-light de¬ 
vice is the bandwidth-delay product, which we estimate 
using the relation rAf = /AA[ng(AT) — ^mim(At)]/At- 
Here, AA is the FWHM of the transmittance peak, 
while the group indices of the metamaterial rig and the 
MIM waveguide timim are evaluated at the wavelength 
of maximum transmittance At- For a single metaatom 
we find the bandwidth-delay product to be 0.26, which 
is quite large considering the short propagation length. 
Joining additional metaatoms significantly increases the 
bandwidth-delay product to 0.36, 0.68, and 1.2 for 2, 4, 
and 8 cells, respectively. Despite the decrease in FWHM 
of the transmittance peak, the increase in propagation 
length and in group index results in an overall increase 
in the bandwidth-delay product. The large values for 
the bandwidth-delay product come at the expense of de¬ 
creasing maximum transmittance. However, even for 8 
cells with a huge bandwidth-delay product of 1.2, the 
maximum transmittance is still 10%. Adding more cells 
increases the bandwidth-delay product even further, but 
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the transmittance level becomes impractical. Given a 
certain tolerance level for the maximum transmittance, 
say 10%, the maximum number of joined metaatoms 
can actually be estimated from the imaginary part of 
the effective index of a single metaatom [Fig. [3][c)] as 
N = — ln(T)/[2Im(neff)a/co]. This is possible due to the 
fast convergence of the effective index in the FIT window 
with increasing number of cells. 

While we only showed the long-wavelength EIT win¬ 
dow in Fig. [3l we have also performed simulations on the 
short-wavelength EIT window of Fig. [2l Using the same 
approach for creating the metamaterial, we find compa¬ 
rable values for the group index and bandwidth-delay 
product, providing support that our slow-light metama¬ 


terial can be made to operate using different Fabry-Perot 
modes at different wavelengths. In fact, the detuned res¬ 
onators of the individual metaatoms can be achieved with 
other Fabry-Perot resonators than the aperture-coupled 
stub-shaped considered here, such as nanodisk cavities^^ 
or T-shaped resonators^. 

In summary, we have considered a simple configura¬ 
tion for designing slow-light metamaterials with large 
bandwidth-delay products based on metaatoms exhibit¬ 
ing EIT. The EIT effect is achieved using detuned res¬ 
onators to mimic the atomic dressed-state picture. By 
adjusting the detuning, the metamaterial can be tuned 
to different operating wavelengths. 
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